The Astrophysical Journal Letters, 717: 

Preprint typeset using I4TeX style emulateapj v. 1 1/10/09 



L118-L121, 2010 July 10 



DOl:10.1088/2041-8205/717/2/LllS 



O 

(N 

1—5 

(N 

o' 

iHj 

>■ 
(N ■ 

o\ 

rn 

vn 

o, 
o: 



% 



GEV BREAKS IN BLAZARS AS A RESULT OF GAMMA-RAY ABSORPTION WITHIN THE BROAD-LINE REGION 

JuRi Poutanen' and Boris Stern'-'-'' 

^Astronomy Division, Department of Physics, P.O.Box 3000, 90014 University of Oulu, Finland; juri.poutanen@oulu.fi 
^Institute for Nuclear Research. Russian Academy of Sciences, Prospekt 60-letiya Oktyabrya 7a, Moscow 1 17312, Russia; boris.stem@gmail.com 
^Astro Space Center, Lebedev Physical Institute, Profsoyuznaya 84/32, Moscow 117997, Russia 
Received 2010 May 20; accepted 2010 June 8; published 2010 June 21 

ABSTRACT 

Spectra of the brightest blazars detected by the Fermi Gamma-ray Space Telescope Large Area Telescope 
cannot be described by a simple power law model. A much better description is obtained with a broken power 
law, with the break energies of a few GeV. We show here that the sharpness and the position of the breaks can 
be well reproduced by absorption of 7-rays via photon-photon pair production on He II Lyman recombination 
continuum and lines. This implies that the blazar zone lies inside the region of the highest ionization of the 
broad-line region (BLR) within a light-year from a super-massive black hole. The observations of 7-ray spectral 
breaks open a way of studying the BLR photon field in the extreme-UV/soft X-rays, which are otherwise hidden 
from our view. 

Subject headings: BL Lacertae objects: general - galaxies: active - galaxies: jets - gamma rays: general - 
radiation mechanisms: non-thermal 



1. INTRODUCTION 

The Fermi Gamma-ray Space Telescope has detected more 
than a hundred of blazars in the 100 MeV-100 GeV range 
JAbdo et alj2009bl) and allowed their detailed spectral studies 
with unprecedented accuracy. A surprising result is that the 
spectra of high-luminosity sources, flat-spectra radio quasars 
(FSRQs) and low-energy synchrotron peaked BL Lac objects, 
are much better described by a broken power la w than by a 
simple powe rlaw or any smoothly curved models dAbdo et alj 
I2009all20r0.) . The break energies, mostly lying in the 2-10 
GeV range (as measured in the object frame), seem to be too 
small to be produced by the 7-ray absorption due to photon- 
photon pair production in the broad-line region (BLR), as 
the strongest BLR line, Lya, absorbs 7-rays starting from 
only 25.6 GeV. The observed spectral breaks are also too 
sharp and the change in the spectral indices are too large to 
be associat ed with the c ooling or the Klein-Nishina effects 
dGhisenini & Tavecchidl2 "009). 

Calculations of the 7-ray absorption within the BLR are 
usually limited to the contribution from lines in the observed 
optical ban d (e.g. [Li u & Bai 2006; Reimer 2007) with a few 
exceptions ( iTavecchio & Mazin.2009.) . Theoretical models, 
however, predict strong helium UV lines a s well as metal 
lines in the soft X-rays (see, e.g. lKroliMll999l) . because a typ- 
ical quasar spec trum has a p ower law tail extending to the 
X-ray band (Laoret a"Ll ll997l) . In addition to the lines, there 
are strong recombination continua of hydrogen and He II at 
13.6 and 54.4 eV, respectively. These sharp, line-like features 
cause jumps in the 7-ray opacity at '^19.2 and 4.8 GeV. We 
propose here that the observed spectral breaks in blazars are 
produced by absorption on photons of these recombination 
continua. 



2. GAMMA-RAY ABSORPTION WITHIN THE BROAD-LINE REGION 

2.1. Emission From Broad-Line Region 

The BLR around quasars is known to emit a number of 
strong lines associated with very different ionization stages. 
Reverberation mapping of Seyfert galaxies and quasars using 
the C IV A 1549 line showed that the size of the BLR scales 



with the luminosity as (iKaspi et alJl2007h 



-Re IV, 18 ~ 0.4L, 



1/2 
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The estimated BLR size is, however, 2-3 times larger for 
the B aimer lines and 3 times smaller for high-ionization 
Hen A 1640 and Nv A1240 lines (Koristaetal. 1991 
iPeterson & Wandelll 19991) . This argues in favor of strong ra- 
dial ionization stratification of the BLR spanning over an or- 
der of magnitude in radius. The BLR spectrum is thus very 
much dependent on the distance from the central ionizing 
source. 

Using spectral synthe sis code XSTAR (version 2.2, 
iKallman & Bautistall2001h . we generated a grid of photoion- 
ization models of BLR clouds assumed to be simple slabs of 
constant gas density and a clear view to the ionizing source. 
For simplicity, we fix the cloud column density at iVn = lO'^'^ 
cm^-^ and vary the ionization parameter ^ = _L/(7-^nH) from 
IqO.'o jq 10^-^. We assume a gradual change of the cloud 
density as a function of distance from the central source, 
nn oc r^^. Taking n^ — 10^^'^ cm^"^ at the inner edge of the 
BLR, the distances r to the ionizing source are between 0.03 
(logi^ — 2.5) and 1 pc (log^ = 0.5) for a quasar luminosity 
L — 10^^ erg s^^ . This BLR model is broadly consistent with 
the co nstraints from reverberation mapping ( Kaspi & Netz^ 
I1999I) . It is, of course, still an oversimplification, but this is 
enough to show the impact of BLR on propagation of the GeV 
photons. 

The quasar spectrum is taken as a sum of the standard multi- 
color accretion disk plus a p ower law of total luminosity 10% 
and extending to 100 keV (iLaor et anil997h . The resulting 
BLR spectra are shown in FigurelHa). Their absolute normal- 
ization depend on the covering fraction of the BLR clouds. If 
it does not vary dramatically over distance, the radiation field 
at a given radius is dominated by the locally produced pho- 
tons corresponding to a specific ionization parameter As our 
aim is not the study of the BLR emission, but its effect on 
the high-energy photon propagation, we smooth the spectrum 
to show the most significant features (see Table [T]i- In the 
high-ionization zone at small distances to the central source. 
He II Ly lines and recombination continuum at 40-60 eV are 
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Fig ure 1. (a) Spectrum of the broad emission line region. The spectra for various ionization parameters § computed using a photoionization code XSTAR version 
2.2 IKallman & Bautista"2001') are averaged on a logarithmic energy grid with the 5% spectral resolution to show the most significant features dominating the 
total flux. The incident typical quasar spectrum (Laor et al. 1997) shown by a smooth curve consists of a multicolor dislc with the maximum temperature 10"" K 
and a power law tail of photon index F = 2 extending to 100 keV with 10% of the total luminosity, (b) Average cross section (in units of Thomson cross section) 
for 7-rays due to photon-photon pair production on the BLR spectra shown in panel (a). Opacity produced by the disk radiation, direct or scattered outside 
the BLR clouds (see e.g. Blandford & Levinson 1995), is neglected here. Jumps in opacity are clearly seen at energies corresponding to the strongest lines and 
recombination continua. The 0.3-0.5 GeV jump for large ^ is produced by O VIII AA16-19 lines, the jump at low g is due to the O VII A22 line complex (see 
Table[T). The jump at ~5 GeV visible at all ionizations is due to the He II recombination continuum and Ly lines at 40-60 eV The 20-30 GeV jump seen at low 
^ is produced by H I Ly lines and continuum and the C IV A 1549 line. 



Table 1 

BLR Strongest Lines and Recombination Continua 
Causing Jumps in the 7-Ray Opacity 



isotropic photons have strong line-Uke features peaking at 
energy Eq, the photon-photon pair production cross sec- 
tion (777(5) for a photon of energy E is the function of 



Feature 



(A) 



(eV) 



(GeV) 



Low-ionization lines 
O VII (blend) 22 560 0.47 

Cv (blend) 40.5 305 0.86 

Low-ionization H I 10 eV complex 
Ly continuum 911 13.6 19.2 

Lyo 1215 10.2 25.6 

Civ 1549 8.0 32.6 

Low-ionization He I 20 eV complex 
He I rec. continuum 504.2 24.6 10.6 

He I 584.3 21.2 12.3 

High ionization lines 
Ovm 16.01 774 0.34 

Ovill 18.97 653 0.40 

Cvi 33.74 367 0.71 

He II 50 eV complex 
He II Ly continuum 227.8 54.4 4.8 

Fexv 284.2 43.6 

Si XI 303.3 40.9 

HellLya 303.8 40.8 6.4 

^ Wavelength of the spectral feature. 

Energy of the spectral feature i?BLR(eV)« 
12393/Ablr(A). 

'^ Energy of the jump in the 7-ray opacity i?^(GeV) fti 
261/Bblr (eV) » Ablr(A)/47.5. 

relatively strong. There are also significant soft X-ray lines 
of O Vlll and C VI, which are important for estimation of the 
^500 MeV photon propagation. In the low-ionization zone at 
larger distances, hydrogen Ly lines and recombination contin- 
uum dominate the photon flux with the additional strong lines 
from C IV, C V, He l. He 11, as well as helium recombination 
continua. 

2.2. Absorption of GeV Photons 

Photons produced in the BLR constitute the target for the 
GeV photons propagating through them. If the background 



the product 01 two energies s = fcfco/ tTTZeC ) (see, e.g. 
iGould & Schrededll967l; IZdziarski|[T988h . It has a thresh- 
old s = 1, rapidly grows to the maximum of about 20% of 
the Thomson cross section ctt at s « 3, and then slowly de- 
creases. The optical depth for a photon of energy E through 
the region of size R filled with photons of column density A^ph 
is Ty'y{E, Eq) — TTajj{s)/aT, where 



tt ~ Npi^cFT = 



L(Tt^ 



= 110- 



L45 10 eV 



(2) 



AnRcEo Ris Eq 

and L is the line luminosity (we defined Q = lO^Q^; in cgs 
units). Equation (|2]i implies that it is not the luminosity, but 
its ratio to the typical size (i.e. compactness) that determines 
a role of a specific line in absorption of the 7-rays. Because of 
the luminosity dependence of the BLR size, the opacity also 
depends on it as 

TT ex L^/^ (3) 

This provides a simple explanation of why the GeV breaks are 
observed only in the brightest blazars. 

The spectrum transmitted through the region is attenuated 
as oc exp{—Tj^{E, Eq)). If the incident spectrum is a power 
law, a break appears above the threshold energy {rricc'^)'^ / Eq, 



where the spectral index changes by 

dlnexp{~T^j{E, Eq)) 



AT 



,(s) 



dln£; 



Ins 



H 
4 



(4) 



The efficiency of 7-ray absorption within a BLR depends 
on the photon distribution, with the averaged photon-photon 
pair production cross section being 

CT-wlii) = ^^— I (J^jis)Nph{Eo)dEo. (5) 



'jj{E) - —— 

JVph 



Figure [Tfb) presents the photon-photon pair production cross 
section averaged over the BLR spectra shown in Figure [Tf a). 



GeV breaks in blazars as a result of gamma-ray absorption within BLR 



Table 2 

Spectral Properties of Blazars 



Object 



Power Law 



Broken Power Law 

r2 Sbrcak(l + z)(GeV) 



Power Law + Double Absorber 



T-Ra 



T"H 



3C 454.3 


0.859 


117 


2.36±0.02 


3.60±0.22 


4.5±0.5 


6.5 


2.37±0.02 


6.1 ±0.9 


18.5+^8 


4.1 


PKS 1502+106 


1.839 


55 


2.15±0.03 


2.87±0.16 


7.8±1.5 


7.8 


2.13±0.03 


1.6±0.6 


8.4±1.6 


6.3 


3C279 


0.536 


18 


2.17±0.07 


2.56±0.09 


1.8±0.6 


4.6 


2.28±0.04 


2.0±1.1 


4.5±3.1 


10.1 


PKS 1510-08 


0.36 


13 


2.43±0.05 


2.84±0.27 


3.1±1.8 


6.6 


2.45±0.04 


2.7±1.5 


2 7+8 
Z./_2 7 


8.1 


3C273 


0.158 


10 


2.82±0.06 


3.40±0.42 


1 o+i.o 

^■' 1.9 


6.1 


2.87±0.05 


3.61L 


0+°° 
" 


7.8 


PKS 0454-234 


1.003 


50 


2.04±0.05 


2.81±0.17 


5.3±1.0 


12.3 


2.04±0.04 


3.0±0.8 


9.5±2.7 


13.7 


PKS 2022-07 


1.388 


15 


2.45±0.05 


3.02±0.17 


9.6±4.3 


11.6 


2.48±0.06 


0.81- 


2.9tt- 


12.9 


TXS 1520+319 


1.487 


11 


2.49±0.07 


2.89±0.24 


4.7±0.5 


7.9 


2.48±0.74 


1.7±1.6 


6.5l« 


7.2 


RGB J0920+446 


2.19 


21 


1.99±0.08 


3.47±0.4 


19±5 


7.8 


2.01±0.07 


0+0.5 

"-0 


7.6±2.9 


11.9 


Note. — The number of degrees 


of freedom is 


12 for the power law model and 10 for other models. 













If the ionization is high, the dominant photon source is the 
He II complex at 40-60 eV (see Table [T]i, which produces a 
break in opacity at 4-7 GeV. Above 10 GeV, the opacity is a 
rather smooth function of energy with additional absorption 
coming from He I, hydrogen Ly lines and recombination con- 
tinuum, and C V A2274. In the 0.3-0.7 GeV region the opacity 
has another break due to high-ionization lines of O VII. De- 
tection of this break requires tt > 300, which certainly is 
enough to completely absorb the GeV radiation. 

In the low-ionization environment, the strongest absorption 
is produced by hydrogen Ly lines and recombination contin- 
uum, forcing a break in opacity at about 20-30 GeV. At the 
same time, the break due to helium is also clearly visible at a 
few GeV. The breaks below 1 GeV would be probably impos- 
sible to detect due to low opacit y at these energies. The opac- 
ity is nearly flat up to 1 TeV (jTavecchio & Mazinll200 9) be- 
cause of the contribution from additional lines: Mg II A2800, 
Ha, andHelA10832. 

The breaks at 4-7 and 20-30 GeV are expected for a large 
range of ionization parameters. The total opacity in 5-50 GeV 
range can therefore be very roughly represented as a sum of 
opacities resulting from two "lines" at ionization energies of 
He II and H I, 54.4 and 13.6 eV, taken in different proportions. 
The opacity in these lines can be characterized by the optical 
depths tho and th, which reflect the column density of pho- 
tons at respective energies via Equation (|2]i. In this simple 
approximation, which we will call a double-absorber model, 
the breaks appear at 4.8 and 19.2 GeV. 

3. BLAZARS OBSERVED BY FERMI 

We analyze the data for several brightest FSRQs from the 
sample of 12 objects in Table 1 in Abdo et al. (2010) and 
choose the same 180 day interval for easier comparison. We 
cannot use the standard software distributed by the LAT team, 
because it does not contain the required spectral models. In- 
stead, we wrote our own version and check its performance, 
by comparing its results for t he simple b r oken power law 
models with those obtained by lAbdo et aTl 42010) using the 
standard maximum likelihood analysis tool gtlike. We ex- 
cluded three objects for which our simplified analysis is not 
adequate: 4C+38.41 and PKS 0528+134, because of severe 
source confusion, and PKS 1908-201 since it is close to the 
Galactic bulge and the background is high. 

3.1. Data Analysis 

We use class 3 (diffuse) photons and impose the cuts at 
zenith angle <105° and incident angle to the detector axis 



<60°. The latter cut excludes counts with very low value 
of the detector response and restricts variation of the count 
weight by factor ~2.7. We used P6_V3_DIFFUSE version 
of the response function. The background was measured in a 
circle with the radius 6°, separated from the object by ^15° 
and avoiding local sources. We accumulate counts in the cir- 
cle centered at the source location with the energy-dependent 
radius r = max{10°(-B/100MeV)-''^0.6°} (which coiTe- 
sponds to 95% containment). 

At this stage, we are not interested in the absolute normal- 
ization of the spectrum, but only in its shape, therefore our 
strategy is to sample the response function at photon arrival 
times (rather than integrating it over time). As the energy res- 
olution in the range of interest (1-30 GeV) is good enough, 
we can treat the response matrix as diagonal and just prescribe 
to each photon the weight inversely proportional to the effec- 
tive area for the given energy and incident angle. We restrict 
the spectral fit to energies above 200 MeV We use binned fit- 
ting as it provides us with a very simple and clear indicator of 
the fit quality, x^. For a few bins at high energies, where the 
number of photons is low, we use Poisson likelihood adding 
—2 logP(n, /i) to x^ (here n is the number of counts in the 
bin and /i is the prediction of the model). The number of 
such bins is small and the meaning of x^ is not significantly 
affected. For the minimization we use the standard code MI- 
NUIT from the CERN library. 

3.2. Spectral Fits 

First, we fit the data with a simple power law model. For 
five objects out of nine from our sample, the fits are acceptable 
(X^ < 18 for 12 degrees of freedom). Next, we apply the bro- 
ken power law model with parameters Fi , r2, -Ebroaic and free 
normalizatio n (see Tablel2ll. Th ere is a good agreement with 
the results bv lAbdoet alTMT^ , except for RGB J0920+446 
(where the best- fit parameters given by Abdo et al. (2010) ev- 
idently do not match the spectrum of the object presented in 
the same paper) and for PKS 0454-234 (where our minimiza- 
tion routine probably found a different local minimum). The 
breaks are not statistically significant in PKS 1510-08, 3C 
273, PKS 2022-07, and TXS 1520+319. 

Finally, we have fitted the data with a power law and the 
double-absorber model taking optical depths thc and th, pho- 
ton index F and normalization as free parameters. The "lines" 
were redshifted by the appropriate 1 + z factor. The fits 
give statistically significant detections of absorption with a 
good x^ for 3C 454.3, PKS 1502+106, PKS 0454-234, and 
RGB J0920+446 (see Table |2] and Figure ^. Upper limits 
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Figure 2. Specti'al energy distribution of a few blazars as observed with 
FennifLPd. The best-fit broken power law and a power law with the double- 
absorber models are shown by the dashed and solid lines, respectively. 

on absorption for PKS 1510-08, PKS 2022-07, and TXS 
1520-1-319 are significantly smaller than the absorption opti- 
cal depth for 3C 454.3. There are no significant constraints 
for 3C 273. The quality of the fits with the absorption model 
is about the same as with the broken power law model with 
the same number of parameters. In fact, the absorption model 
is less flexible as the energy, where the power law spectrum 
breaks, is fixed. 

4. DISCUSSION AND SUMMARY 

The GeV breaks observed in blazars are well described 
by 7-ray absorption via photon-photon pair production on 
Hell and Hi recombination continuum photons. In RGB 
J0920-I-446, the absorption is seen only at high energies with 
the break energy -Ebicak corresponding to the pair-production 
threshold on hydrogen recombination photons, while in other 
cases the break is close to the threshold for the absorption 
on He II recombination continuum. The exact position of the 
break depends on the ionization parameter that determines the 
contribution of metals and affects the position of the centroid 
of the 50 e V complex. If thc is small, then the break shifts to 
19 GeV as observed in RGB J0920H-446. 

A rather large ratio of the fitted optical depths thc/th ^ 
1/4, implies that the 7-ray emitting region has to lie within 
the high-ionization zone of the BLR with log^ > 2. For the 
brightest 7-ray object in our sample, 3 C 454.3 , with its accre- 
tion luminosity of about 10''^ erg s~^ (I Smith et al...l98 8). the 
high-ionization zone should be within about 0.1 pc. This cor- 
responds to about lO"^ Schwarzschild r adii for an ~ IQ^Mq 
central black hole jBonnoli et alj|2010l) . At such a distance, 
the luminosity in the 50 eV complex can be as small as '^10'*'* 
erg s^^, a per mille of the accretion luminosity, to provide the 
necessary opacity with thc ^ 6 (see Equation^. 

The opacity measured above 20 GeV and the Lyg luminos- 
ity of IQ^^ erg s^^ observed in 3C 454.3 (IWills et al.lll995h 
allow us to estimate the Lya emission zone size. Taking the 
recombination continuum luminosity equal to that of Lya and 



using Equation dU we get i?LyQ = 2 ± 1 pc, which indicates 
that absorption at these energies happens at a larger distance 
than absorption by He II photons. 

The constraints obtained on the 7-ray emission site im- 
ply that the jet is already accelerated to a relativistic veloc- 
ity within a thousand gravitational radii. It also strongly 
constrains the mechanisms for 7-ray production. The pos- 
sible sources of soft photons for Comptonization i n the jet 
are th e accretion disk (jDermer & Schhckeiserlll993h and the 
BLR (iSikora et al.l ll994). As a source of photons, the dust 
emiss ion at 10 pc scale (Blazejowski et al. 2000; Sikora et al] 
l2008h cannot be important. 

Let us also remark that the GeV photons absorbed in the 
BLR produce electron-positron pairs which, spiraling in the 
magnetic field, radiate away their energy isotropically. This 
emission cannot compete with the beamed emission from the 
jet, but can contribute to the high-energy emission of radio 
galaxies observed at large angles to the jet axis (see, e.g. 
iRoustazadeh Sh. & B oettcher 2010). 

Our interpretation of the GeV breaks implies that additional 
breaks (depressions) at 0.3-0.7 GeV produced by the soft X- 
ray lines within the high-ionization zone should be seen, once 
the photon statistics is high enough. The 7-ray spectroscopy 
can be used as a powerful tool for studying the extreme-UV 
and soft X-ray emission in the quasars' vicinity, which is oth- 
erwise hidden from us by interstellar absorption. 
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